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Abstract 

The DAMA/Nal experiment (~ 100 kg highly radiopure Nal(Tl)) was pro- 
posed, designed and realised to effectively investigate in a model independent 
way the presence of a Dark Matter particle component in the galactic halo by 
exploiting the annual modulation signature. With a total exposure of 107731 
kg • day, collected over seven annual cycles deep underground at the Gran Sasso 
National Laboratory of the I.N.F.N., it has pointed out - at 6.3 o C.L. - an effect 
which satisfies all the peculiarities of the signature and neither systematic effects 
nor side reactions able to mimic the signature were found. Moreover, several (but 
still few with respect to the possibilities) corollary model dependent quests for 
the candidate particle have been carried out. In this paper the obtained results 
are summarized and some perspectives are discussed at some extent. 

Keywords: Dark Matter; WIMPs; underground Physics 
PACS numbers: 95.35. +d 

1 Introduction 

The problem of the existence of Dark Matter in our Universe dates back to the astro- 
physical observations at the beginning of past century [1, 2], but the presence of Dark 
Matter in our Universe has been definitively accepted by the scientific community only 
about 40 years later, when two groups performed systematic measurements of the ro- 
tational velocities of celestial bodies in spiral galaxies [3]. After the '70 many other 
observations have further confirmed the presence of Dark Matter in the Universe and, 
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at present, the measurements are mainly devoted to the investigation of the quan- 
tity, of the distribution (from the cosmological scale down to the galactic one) and of 
the nature of the Dark Matter in the Universe. Recent measurements of the CMB 
temperature anisotropy by WMAP [4], analysed in the framework of the Big Bang 
cosmological scenario, support a density of the Universe: = 1, further crediting that 
most of the Universe is dark. Recently, it has been suggested from observations on 
the supernovae la at high red-shift as standard candles that about 73% of fl might 
be in form of a dark energy [5]; the argument is still under investigation and presents 
some problems on possible theoretical interpretations. However, even in this scenario 
- where the matter density in the Universe would be 0^ ~ 0.3 [4, 6] - large space 
for Dark Matter particles in the Universe exists. In fact, the luminous matter can 
only account for a density ~ 0.005 and the baryonic Dark Matter for ~ 0.04. On the 
other hand, the contribution of Dark Matter particles relativistic at the decoupling 
time is also restricted to be < 0.01 by considerations on large scale structure forma- 
tions [7]. Thus, most of the Dark Matter particles in the Universe, relics from the 
Big Bang, were non relativistic at decoupling time; they are named Cold Dark Matter 
particles (CDM). The CDM candidates have to be neutral, stable or quasi-stable (e.g. 
with a time decay of order of the age of the Universe) and have to weakly interact 
with ordinary matter. These features are respected by the axions (also investigated 
by DAMA/Nal [8]) and by a class of candidates named WIMPs (Weakly Interacting 
Massive Particles). In particular, in the Standard Model of particle Physics no particle 
can be a suitable candidate as CDM; thus, a new window beyond the Standard Model 
can be investigated. 

Finally, several observations have suggested that our Galaxy should also be cm- 
bedded in a dark halo with mass at least 10 times larger than that of the luminous 
matter. The verification of the presence of a Dark Matter particle component in the 
galactic halo has been the main goal of the DAMA/Nal experiment. 

2 The DAMA/Nal experiment 

The DAMA/Nal experiment was proposed in 1990 [9], designed and realized having 
the main aim to investigate in a model independent way the presence of a Dark Matter 
particle component in the galactic halo [10, 11, 12, 13, 14, 15, 16, 17, 18]. For this 
purpose, we planned to exploit the effect of the Earth revolution around the Sun on 
the Dark Matter particles interactions on the target-nuclei of suitable undergroimd 
detectors. In fact, as a consequence of its annual revolution, the Earth should be 
crossed by a larger flux of Dark Matter particles in June (when its rotational velocity 
is summed to the one of the solar system with respect to the Galaxy) and by a smaller 
one in December (when the two velocities are subtracted) (see Fig. 1). This offers an 
efficient model independent signature, able to test a large interval of cross sections and 
of halo densities; it is named annual modulation signature and was originally suggested 
in the middle of '80 by [19]. 

In fact, the expected differential rate as a function of the recoil energy, dR/dEn 
(see rcf. [18] for detailed discussion), depends on the WIMP velocity distribution and 
on the Earth's velocity in the galactic frame, Ve{t)- Projecting Ve{t) on the galactic 
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-30 km/s 



Figure 1: Schematic view of the Earth motion around the Sun. 



plane, one can write: Ve{t) — Vq + Vt^cos^cosuj{t — to). Here ^0 is the Sun's velocity 
with the respect to the galactic halo [vq ~ + 12 km/s and is the local velocity 
whose value is in the range 170-270 km/s [12, 20]); = 30 km/s is the Earth's orbital 
velocity around the Sun on a plane with inclination 7 = 60° with the respect to the 
galactic plane. Furthermore, lj= 27r/T with T=l year and roug hly to ~ 2'"^ June 
(when the Earth's speed is at maximum). The Earth's velocity can be conveniently 
expressed in unit of vq: r]{t) = Ve(t)/vQ — rjo + Ar/cosuj{t — tg), where - depending on 
the assumed value of the local velocity - ?7o = l-04-1.07 is the yearly average of rj and 
Ar] = 0.05-0.09. Since Arj <C tjo, the expected counting rate can be expressed by the 
first order Taylor approximation: 

dR , , dR , , d f dR\ 

Averaging this expression in a /c-th energy interval one obtains: 

drj 



Skivit)] = 'S'fc[77o] + [-^],,oA77cosw(t - to) = So^k + Sm,kCosuj{t - to), (2) 



with the contribution from the highest order terms less than 0.1%. Thus, the annual 
modulation signature is very distinctive since a WIMP-induced seasonal effect must 
simultaneously satisfy all the following requirements: (i) the rate must contain a com- 
ponent modulated according to a cosine function; (ii) with one year period; (iii) a 
phase that peaks roughly around ~ 2"'' June; (iv) this modulation must only be found 
in a well-defined low energy range, where WIMP induced recoils can be present; (v) it 
must apply to those events in which just one detector of many actually "fires", since 
the WIMP multi-scattering probability is negligible; (vi) the modulation amplitude in 
the region of maximal sensitivity must be <7% for usually adopted halo distributions, 
but it can be larger in case of some possible scenarios such as e.g. those in refs. [21, 22]. 
Only systematic effects able to fulfil these 6 requirements and to account for the whole 
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observed modulation amplitude could mimic this signature; thus, no other effect in- 
vestigated so far in the field of rare processes offers a so stringent and unambiguous 
signature. With the present technology, the annual modulation signature remains the 
main signature of a WIMP signal. 

Of course, the amount of the measured effect depends e.g.: i) on the sensitivity 
of the experiment to the coupling and to the particle physics features of the WIMP; 
ii) on the nuclear features of the used target-nucleus; iii) on the features of the dark 
halo and, in particular, on the spatial and velocity distributions of the Dark Matter 
particles. The quality of the running conditions also plays a crucial role to detect such 
a rare process. 

Considering its main goal, DAMA/Nal has been designed by employing and further 

developing all the necessary low-background techniques and procedures. A detailed 
description of the set-up, of its radiopurity, of its performance, of the used hardware 
procedures, of the determination of the experimental quantities and of the data reduc- 
tion has been given in refs. [23, 13, 14, 18]. Here only few arguments are addressed. 

The DAMA/Nal experiment was located deep underground in the Gran Sasso 
National Laboratory of I. N.F.N. ^, whose main features have been reported in [27, 28, 
29, 30]. 

In particular, the Nal(Tl) scintillator was chosen as the best target material to 
investigate the process since it offers e.g.: 

• well known technology 

• reachable high radiopurity by material selections and protocols, by chemi- 
cal/physical purifications, etc. 

• large mass 

• high duty cycle 

• feasible well controlled operational conditions and monitoring 

• routine calibrations feasible down to keV range in the same conditions as the 
production runs 

• high light response, that is keV threshold reachable 

• absence of the necessity of re-purification or cooling down/ warming up proce- 
dures (implying high reproducibility, high stability, etc.) 

• absence of microphonic noise and an effective noise rejection at threshold (time 
decay of Nal(Tl) pulses is hundreds ns, while that of noise pulses is tens ns) 

• sensitivity to spin- independent (SI), spin-dependent (SD) and mixed (SI&SD) 

couplings as well as to several other existing scenarios 

•^Wc take this occasion to remind that DAMA/Nal has been part of the DAMA project, which is 
also coinijosed by several other low background set-ups, such as: i) DAMA/LXe (~ 6.5 kg pure liquid 
Xenon scintillator) [24, 25]; ii) DAMA/R&D, set-up devoted to tests on prototypes and small scale 
experiments [26]; iii) the new second generation large mass Nal(Tl) radiopure set-up DAMA/LIBRA 
(see later); iv) DAMA/Ge detector for sample measurements. 
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• sensitivity to both high (by Iodine target) and low (by Na target) mass candidates 

• possibility to effectively investigate the annual modulation signature in all the 
needed aspects 

• pulse shape discrimination feasible at reasonable level 

• possibility to achieve significant results on several other rare processes 

• no safety problems 

• necessity of a relatively small underground space 

• the lowest cost with the respect to every other considered technique 

However, neither commercial low background Nal(Tl) detectors nor Nal(Tl) detec- 
tors grown with old technology (even after " revision" ) can reach the needed sensitivity; 
thus, a devoted special R&D was realized with the Crismatec company. In this frame- 
work, all the needed materials were selected and through devoted and severe protocols 
the final DAMA/Nal highly radiopure detectors were realized ^. 

The performances of the nine 9.7 kg highly radiopure DAMA/Nal detectors, in- 
cluding their measured radiopurity, arc discussed in details in [23, 14, 18]. The bare 
Nal(Tl) crystals are encapsulated in radiopure Cu housings; moreover, 10 cm long 
Tetrasil-B light guides act as optical windows on the two end faces of each crystals 
and are coupled to specially developed low background photomultipliers (PMT). The 
measured light response is 5.5 - 7.5 photoelectrons/keV depending on the detector. 
The two PMTs of a detector work in coincidence with hardware thresholds at the sin- 
gle photoelectron level in order to assure high efficiency for the coincidence at few keV 
level. The energy threshold of the experiment, 2 keV, is instead determined by means 
of X-rays sources and of keV range Compton electrons on the basis also of the features 
of the noise rejection procedures and of the efficiencies when lowering the number of 
available photoelectrons [23]. 

The detectors are enclosed in a sealed copper box, continuously maintained in HP 
Nitrogen atmosphere in slightly overpressure with respect to the external environment. 

A suitable low background hard shield against electromagnetic and neutron back- 
ground was realized using very high radiopure Cu and Pb bricks [23], Cd foils and 
10/40 cm polyethylene/paraffin; the hard shield is also sealed in a plexiglas box and 
maintained in the high purity (HP) Nitrogen atmosphere. Moreover, about 1 m con- 
crete (made from the Gran Sasso rock material) almost fully surrounds the hard shield 
outside the barrack and at its bottom, acting as a further neutron moderator. 

A three-level sealing system from environmental Radon is effective. In fact, the 
inner part of the barrack, where the set-up is allocated, has the floor (above the con- 
crete) and all the walls sealed by Supronyl (permeability: 2 • 10^^^ cm^/s [31]) plastic 
and the entrance door is air-tight. A low level oxygen alarm informs the operator 
before entering the inner part of the barrack since the HP Nitrogen which fills both 

^Wc remark that the low background technique requires very long and accurate work for the 
selection of low radioactive materials by sample measurements with HP-Ge detectors (placed deep 
underground in suitable hard shields) aud/or by mass spectrometer analyses. Thus, these measure- 
ments are often difficult experiments themselves, depending on the required level of radiopurity. 
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the inner Cu box and the external plexiglas box is released in this closed environment. 
The Radon level inside the (sealed) barrack is continuously monitored and recorded 
with the production data [10, 11, 23, 13, 14, 18]. 

On the top of the shield a glove-box (also maintained in the HP Nitrogen atmo- 
sphere) is directly connected to the inner Cu box, housing the detectors, through 
Cu pipes. The pipes arc filled with low radioactivity Cu bars (covered by 10 cm of 
low radioactive Cu and 15 cm of low radioactive Pb) which can be removed to allow 
the insertion of radioactive sources for calibrating the detectors in the same running 
condition, without any contact with external environment [23]. 

The whole installation is air-conditioned and the operating temperature as well 
as many other parameters are continuously monitored and acquired with the produc- 
tion data. Moreover, self-controlled computer processes automatically monitor several 
parameters and manage alarms [23, 14, 18]. 

The electronic chain and the data acquisition system operative up to summer 2000 
have been described in rcf. [23] , while the new electronics and DAQ installed in summer 
2000 have been described in ref. [18]. 

As regards other aspects, we recall that the linearity and the energy resolution 
of the detectors have been investigated using several radioactive sources [23, 14] such 
as, for the low energy region, ^^Fe (5.9 keV X-rays), ^°^Cd (22 keV X-rays and 88 
keV 7 line) and ^^^Am (59.5 keV 7 line) sources. In particular, in the production 
runs, the knowledge of the energy scale is assured by periodical calibrations with 
^''^Am source and by monitoring (in the production data themselves summed every 
~ 7 days) the position and energy resolution of the 46.5 keV 7 line of the ^^°Pb 
[23, 32, 10, 11, 13, 14, 18]. The latter peak was present at level of few counts 
per day per kg (cpd/kg) - in the measured energy distributions mainly because of a 
contamination (by environmental Radon) of the external surface of the crystals' Cu 
housings, occurred during the first period of the imderground storage of the detectors. 

The only procedure applied to the data regards the noise rejection, which is par- 
ticularly efficient since the difference in time decay of the Nal(Tl) scintillation pulses 
and of the noise pulses is of order of hundreds ns; the procedure have been described 
e.g. in refs. [23, 13, 14]^. 

All the periodical long calibration procedures [23, 32] and the time specifically 
allocated for maintenance and/or for improvements are the main components affecting 
the duty cycle of the experiment. Moreover, in the first period data have been taken 
only in the two extreme conditions for the annual modulation signature. 

The energy threshold, the PMT gain and the electronic line stability are contin- 
uously verified and monitored during the data taking by the routine calibrations, by 
the position and energy resolution of the ^^°Pb line (see above) and by the study of 
the hardware rate behaviours with time. 

Finally, we remind that the experiment took data up to MeV energy region despite 
the optimization was done for the keV energy range. 

The measured low energy distributions of interest for the WIMP investigation have 

®This procedure assures also the rejection of any possible contribution either from afterglows (when 
not already excluded by the dedicated 500 /iS veto time; see above) induced by high energy events or 
from any possible Cerenkov pulse in the light guide or in the PMTs; in fact, they also have fast time 
decay (of order of tens ns) . 
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been given in refs. [33, 13, 14, 34], where the corrections for efficiencies and acquisition 
dead time have already been applied. We note that usually in DAMA/Nal the low 
energy distributions refer to those events where only one detector of many actually 
fires (that is, each detector has all the others in the same installation as veto; this 
assures a further background reduction, which is of course impossible when a single 
detector is used), single-hit events. 

Moreover, the data taking of each annual cycle has been started before the expected 
minimum of the signal rate (which is roughly around ~ 2""^ December) and concluded 
after the expected maximum (which is roughly around ~ 2"'' June). 

The DAMA/Nal set-up has exploited the WIMP annual modulation signature over 
seven annual cycles [10, 11, 12, 13, 14, 15, 16, 17, 18] and the following part of this 
paper will summarize the final model independent result, some of the corollary quests 
for the candidate particle, implications and pespectives. However, it is worth to remind 
that - thanks to its radiopurity and features - DAMA/Nal has also investigated other 
approaches for WIMPs in ref. [32, 33] and several other rare processes such as: possible 
processes violating the Pauli exclusion principle [35], CNC processes in ^^Na and ^^^1 
[36], electron stability [34], searches for neutral SIMPs [37], for neutral nuclearites [37], 
for Q-balls [38] and for solar axions [8]; moreover, DAMA/Nal allowed also the study 
of nuclear transition to a possible superdense nuclear state [39] and possible two-body 
cluster decay of ^^'^I [40] . 

3 The 6.3 a C.L. model-independent evidence for a 
WIMP component in the galactic halo 

As reported in ref. [18], a model independent investigation of the annual modulation 
signature has been realized by exploiting the time behaviour of the residual rates of 
the single-hit events in the lowest energy regions over the seven annual cycles (total 
exposure: 107731 kg • day) [10, 11, 12, 13, 14, 15, 16, 17, 18]. These experimental 
residual rates of the single-hit events are given by: < Vijk — flatjk >jk, where r^fc is 
the rate in the considered i-th time interval for the j-th detector in the fc-th considered 
energy bin and flatjk is the rate of the j-th detector in the fc-th energy bin averaged 
over the cycles. The average is made on all the detectors {j index) and on all the 
energy bins in the considered energy interval. 

This model independent approach offers an immediate evidence of the presence of 
an annual modulation of the rate of the single-hit evcints in the lowest energy region 
over the seven annual cycles as shown in Fig. 2, where the energy and time behaviours 
of the single-hit residual rates are depicted. 

In particular, the data favour the presence of a modulated cosine-like behaviour at 
6.3 a C.L.; in fact, their fit for the (2-6) keV larger statistics energy interval offers a 
modulation amplitude equal to (0.0200 ±0.0032) cpd/kg/keV, a phase: to = (140 ±22) 
days and a period: T = (1.00 ± 0.01) year, all parameters kept free in the fit. The 
period and phase agree with those expected in the case of a WIMP induced effect (T 
= 1 year and to roughly at ~ 152.5-th day of the year). The test on the (2-6) keV 
residual rate in Fig. 2 disfavours the hypothesis of unmodulated behaviour giving a 
probability of 7- 10~^ {x^/d.o.f. = 71/37). Moreover, if the experimental residuals of 
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Fig. 2 are fitted fixing the period at 1 year and the phase at 2"'' June, the fohowing 
modulation ampHtudes are obtained: (0.0233 ±0.0047) cpd/kg/keV, (0.0210 ±0.0038) 
cpd/kg/keV and (0.0192 ± 0.0031) cpd/kg/keV for the (2-4), (2-5) and (2-6) keV 
energy intervals, respectively. 

In Fig. 3 the single-hit residual rate in a single annual cycle from the total exposure 
of 107731 kg • day is presented for two different energy intervals; as it can be seen the 
modulation is clearly present in the (2-6) keV energy region, while it is absent just 
above. The same conclusion is obtained by investigating the data by means of the 
Fourier analysis (performed according to ref. [41] including also the treatment of the 
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Figure 2: Experimental residual rate for single-hit events, in the (2-4), (2-5) and (2-6) 
keV energy intervals as a function of the time elapsed since January 1-st of the first 
year of data taking. The experimental points present the errors as vertical bars and 
the associated time bin width as horizontal bars. The superimposed curves represent 
the cosinusoidal functions behaviours expected for a WIMP signal with a period equal 
to 1 year and phase at 2""^ June; the modulation amplitudes have been obtained by 
best fit. See text. The total exposure is 107731 kg • day. 
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experimental errors and of the time binning); in fact, the results shown in Fig. 4 show 
a clear peak for a period of 1 year in the (2-6) keV energy interval, while it is absent 
in the energy interval just above. 
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Figure 3: Experimental single-hit residual rate in a single annual cycle from the total 
exposure of 107731 kg • day for the (2-6) keV and (6-14) keV energy intervals. The 
experimental points present the errors as vertical bars and the associated time bin 
width as horizontal bars. The initial time is taken at August 7*'*. Fitting the data 
with a cosinusoidal function with period of 1 year and phase at 152.5 days, the follow- 
ing amplitudes are obtained: (0.0195 ± 0.0031) cpd/kg/keV and -(0.0009 ± 0.0019) 
cpd/kg/keV, respectively. Thus, a clear modulation is present in the lowest energy 
region, while it is absent just above. 



A quantitative investigation of the whole energy spectrum up to MeV energy region 
has not shown modulation in any other energy interval (see e.g. [13, 14, 18] and few 
arguments given later). 

Finally, in order to show if the modulation amplitudes are statistically well dis- 
tributed in all the crystals, in all the annual cycles and in the energy bins, the distri- 
butions of the variable S,t^-<S^> reported in Fig. 5. The Sm arc the experimental 
modulation amplitudes for each detector, for each annual cycle and for each considered 
energy bin (taken as an example equal to 0.25 keV), a are their errors and the < Sm > 
represent the mean values of the modulation amplitudes over the detectors and the 
annual cycles for each energy bin. The left panel of Fig. 5 shows the distribution 
referred to the region of interest for the observed modulation, (2-6) keV, while the 
right panel includes also the energy region just above, (2 14) kcV. Since this variable 
is distributed as a gaussian with an unitary standard deviation, the modulation ampli- 
tudes are statistically well distributed in all the crystals, in all the data taking periods 
and in all consiclercxi cnc^rgy bins. 

In conclusion, the data satisfy all the peculiar requirements, given above, for the 
WIMP model independent annual modulation signature. As mentioned, to mimic this 
signature, systematic effects or side reactions shoiild not only be able to fully account 
for the measured modulation amplitude, but also to satisfy the peculiar requirements 
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Figure 4: Power spectrum of the measured single-hit residuals for the (2-6) keV (con- 
tinuous Hne) and (6-14) keV (dotted line) energy intervals calculated according to ref. 
[41] , including also the treatment of the experimental errors and of the time binning. 
As it can be seen, the principal mode present in the (2-6) keV energy interval corre- 
sponds to a frequency of 2.737 • 10"'^ d~^, that is to a period of ~ 1 year. A similar 
peak is not present in the (6-14) keV energy interval. 
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Figure 5: Distributions of the variable 
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the Sm) evaluated for each detector, for each annual cycle and each considered energy 
bin: i) in the region of interest for the observed modulation, (2-6) keV (left panel); ii) 
including also the energy region just above, (2-14) keV (right panel). See text. 



as for a WIMP induced effect. A careful investigation of all the known possible sources 
of systematics and side reactions has been regularly carried out by DAMA/Nal and 
published at time of each data release [10, 11, 13, 14, 18]. 

No effect able to mimic the signature has been found; thus, upper limits (90% 
C.L.) on the possible contributions to the modulated amplitude have been calculated 
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as summarized in Table 1; for a detailed quantitative discussion see ref. [14, 18]. 
In particular, they cannot account for the measured modulation amplitude because 

Tabic 1: Summary of the results obtained by investigating the possible sources of 
systematics or of side reactions [18]. No systematics or side reaction has been found 
able to give a modulation amplitude different from zero; thus cautious upper limits 
(90% C.L.) have been calculated and are shown here in terms of the mcasmcd model 
independent modulation amplitude, S!^^ . As it can be seen, none of them nor their 
cumulative effect could account for the measured modulation; moreover, as discussed 
in details in ref. [14, 18], they cannot mimic the signature. 



Source 


Cautious upper limit 
(90%C.L.) 


Radon 


< 0.2%^^^^ 


Temperature 


< O.bVoST 


Noise 


< l^'S^"' 


Energy scale 




Efficiencies 


< i%sT 


Background 


< Q.h%ST 


Side reactions 


< 0.3%ST 


In addition; no effect can mimic the signature 



quantitatively not relevant and unable to satisfy all the peculiar requirements of the 
signature. 

For the sake of completeness, we remind that possible diurnal effects - correlated 
both with the sidereal and with the solar time - have already been excluded by the 
analysis reported in [33] . 

In particular, as mentioned above, no modulation has been observed in the back- 
ground; in fact, the whole energy spectrum up to MeV energy region has been analyzed 
and the presence of a background modulation in the whole energy spectrum has been 
excluded at a level much lower than the effect found in the lowest energy region [14, 18]. 
This result already accounts also for the background component due to the neutrons 
and the Radon; nevertheless, further additional independent and cautious analyses to 
estimate their possible contribution have been given in ref. [14, 18]. In fact, it has 
been demonstrated in [14, 18] that a modulation of neutron flux - possibly observed 
by the ICARUS coll. as reported in the ICARUS internal report TM03-01 - cannot 
quantitatively contribute to the DAMA/Nal observed modulation amplitude, even if 
the neutron flux would be assumed to be 100 times larger than measured at LNGS 
by several authors with different techniques over more than 15 years. Moreover, in no 
case the neutrons can mimic the signature since some of the peculiar requirements of 
the signature would fail. Similarly, any possible effect of the muon flux modulation re- 
ported by the MACRO experiment [27] is excluded both by quantitative investigation 
[14, 18] and by inability to fulfil all the peculiarities of the signature. 

As regards the possibility of a contribution from the Radon, we remind that the 
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DAMA/Nal has three levels of insulation from the environmental air (see above). 
Moreover, the Radonmeter which continuosly recorded the Radon level inside the 
barrack (that is external to the hard shield and to the detectors) typically measured 
values at level of its sensitivity and no modulation has been observed [10, 11, 13, 14, 
18]. To be on the safest side, even the possible presence of Radon trace in the HP 
Nitrogen atmosphere inside the Cu box, housing the detectors, has been investigated 
by searching for the double coincidences of the gamma-rays (609 and 1120 keV) from 
^^^Bi Radon daughter, obtaining an upper limit of: < 4.5 • 10~^ Bq/m^ (90% C.L.). 
It gives rise to the upper limit reported in Table 1 when assuming an hypothetical 
10%, modulation of possible Radon in the HP Nitrogen atmosphere of the Cu box. 
Anyhow, it is worth to remark that in every case even a sizeable quantity of Radon 
nearby a detector cannot mimic the WIMP annual modulation signature since some 
of the peculiarities of the signature would fail. 
For more detailed discussion see ref. [14, 18]. 

To perform a further relevant investigation, in 1999 we proposed to renew the 
electronic chain of DAMA/Nal removing the multiplexer system and equipping each 
detector with its own transient digitizer. This occurred in summer 2000, thus in the 
last two annual cycles the multiple-hits events and the single-hit events have been 
acquired and analyzed using the same identical hardware and the same identical soft- 
ware procedures. The multiple-hits events class - on the contrary of the single-hit 



Figure 6: Experimental residual rates over seven annual cycles for single-hit events 
(open circles) - class of events to which WIMP events belong - and over the last two 
annual cycles for multiple-hits events (filled triangles) - class of events to which WIMP 
events do not belong - in the (2-6) keV cumulative energy interval. They have been 
obtained by considering for each class of events the data as collected in a single annual 
cycle and using in both cases the same identical hardware and the same identical 
software procedures. The initial time is taken on August 7*''. See text. 

one - does not include events induced by WIMPs since the probability that a WIMP 
scatters off more than one detector is negligible. The obtained results are depicted 
in Fig. 6. The fitted modulation amplitudes arc: A = (0.0195 ± 0.0031) cpd/kg/kcV 
and A = — (3.9± 7.9) • 10~^ cpd/kg/keV for single-hit and multiple-hits residual rates. 
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respectively. Thus, evidence of annual modulation is present in the single-hit resid- 
uals (events class to which the WIMP-induced recoils belong), while it is absent in 
the multiple-hits residual rate (event class to which only background events belong). 
Since the same identical hardware and the same identical software procedures have 
been used to analyse the two classes of events, the obtained result offers an additional 
strong support for the presence of a Dark Matter particle component in the galactic 
halo further excluding any side effect either from hardware or from software procedures 
or from background. 

In conclusion, the presence of a Dark Matter particle component in the galactic halo 
is supported by DAMA/Nal at 6.3 a C.L. and the modulation amplitude measured 
over the 7 annual cycles in Nal(Tl) at the location of the Gran Sasso Laboratory for 
the (2 - 6) keV energy region is (0.0200±0.0032) cpd/kg/keV. This is the experimental 
result of DAMA/Nal. It is model independent; no other experiment whose result can 
be directly compared with this one is available so far in the field of Dark Matter 
investigation. 

4 Corollary results: quests for a candidate particle 
in some of the possible model frameworks 

On the basis of the obtained model independent result, corollary investigations can 
also be pursued on the nature and coupling of the WIMP candidate. This latter 
investigation is instead model dependent and - considering the large uncertainties 
which exist on the astrophysical, nuclear and particle physics assumptions and on 
the parameters needed in the calculations - has no general meaning (as it is also 
the case of exclusion plots and of the WIMP parameters evaluated in the indirect 
detection experiments). Thus, it should be handled in the most general way as we have 
preliminarily pointed out with time passing in the past [10, 11, 12, 13, 14, 15, 16, 17] 
and we have discussed in some specific details in ref. [18]. 

It is worth to note that the results presented in the following are, of course, not 
exhaustive of the many possible scenarios which at present level of knowledge cannot 
be disentangled. Some of the open questions are: i) which is the right nature for the 
WIMP particle; ii) which is its right couplings with ordinary matter (mixed SI&SD, 
purely SI, purely SD or preferred inelastic); iii) which are the right form factors and 
related parameters for each target nucleus; iv) which is the right spin factor for each 
target nucleus: v) which are the right scaling laws (let us remind that even for the 
neutralino case in a MSSM framework with purely SI interaction the scenario could 
be drastically modified as pointed out in ref. [42]); vi) which is the right halo model 
and related parameters; vii) which are the right values of the experimental parameters 
within their uncertainties; etc. 

As regards, in particular, the Dark Matter particle-nucleus elastic scattering, the 
differential energy distribution of the recoil nuclei can be calculated [32, 43, 18] by 
means of the differential cross section of the WIMP-nucleus elastic processes given by 
the sum of the SI and the SD contributions. In the SI case, the nuclear parameters can 
be decoupled from the particle parameters and the nuclear cross sections are usually 
scaled to a defined point-like SI Dark Matter particle- nucleon cross section, erg/. In 
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the SD case the notations [15]; tgO = can be used, where the effective 

WIMP-nucleon coupHng strengths for SD interactions. The mixing angle 9 is defined 
in the [0, tt) interval; in particular, 6 values in the second sector account for Op and 
a„ with different signs. Also in the SD case all the nuclear cross sections are usually 
scaled to a defined point-like SD Dark Matter particle- nucleon cross section, asD [18]. 
Thus, the energy distribution of the recoil rate can be written as a function of (Jsi, 
asD and 6. Another important parameter is the local density of the Dark Matter 
particles, pw = ^Po, where po is the local halo density and ^ < 1) is the fractional 
amount of local density of Dark Matter particles. 

Among the ingredients entering in the model dependent analyses there are the 
nuclear SI and SD form factors, which generally depend on the nature of the involved 
particle. In ref. [18] the existing uncertainties in the usually adopted formulation 
for the SI case have been discussed as well as those for the SD one. It has also 
shown that the SD case is even more uncertain since the nuclear and particle physics 
degrees of freedom cannot be decoupled and a dependence on the assumed nuclear 
potential exists. Moreover, further significant uncertainties in the evaluation of the SD 
interaction rate also arise from the adopted spin factor for the single target-nucleus. 
In fact, the available calculated values are well different in different models and, in 
addition, at fixed model they depend on 6. Thus, not only the target nuclei should 
have spin different from zero (for example, this is not the case of Ar isotopes) to be 
sensitive to Dark Matter particles with a SD component in the coupling, but also well 
different sensitivities can be expected among odd-nuclei having an unpaired proton 
(as e.g. ^^Na and ^^''I) and odd-nuclei having an unpaired neutron (as e.g. the odd 
Xe and Te isotopes and ^"^Ge). 

Another scenario also considered in the corollary DAMA/Nal model dependent 
analyses is that of WIMPs with preferred inelastic scattering which has been suggested 
by [21]. In this case the Dark Matter particles can only inelastically scatter off nuclei 
going to excited levels with a S mass splitting. A specific model featuring a real 
component of the sneutrino (for which the mass splitting naturally arises) has been 
given in ref. [21]. It has been shown that for this inelastic scattering a kinematical 
constraint exists which favours target-detectors media with heavy nuclei (such as ^^^I) 
with the respect to those with lighter ones (such as e.g. ""*Ge). In fact, this process 
can only occur if the particle velocity is larger than a threshold value; this kinematical 
constraint becomes increasingly severe as the nucleus mass decreases [21]. Moreover, 
this model scenario implies some interesting peculiar features when exploiting the 
annual modulation signature; in fact - with the respect to the case of WIMP elastically 
scattering - it would give rise to an enhanced modulated component of the signal with 
respect to the unmodulated one and to largely different behaviours with energy for 
both the components (showing both higher mean values) [21]. The differential energy 
distribution of the recoil nuclei in the case of inelastic processes is function of ^ap, mw 
and S, analogoulsy as above other ingredients, as e.g. the form factor, play also a role 
[18]. 

As mentioned, the expected energy distribution for the scatterings of Dark Matter 

particles depends on pw and on the velocity distribution of the Dark Matter particles at 
Earth's position. The experimental observations regarding the dark halo of our Galaxy 
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do not allow to get information on them without introducing a model for the Galaxy 
matter density. The dark halo model widely used in the calculations carried out in 
the field of direct detection approaches is the naive isothermal sphere that corresponds 
to a spherical infinite system with a flat rotational curve. Despite its simplicity has 
favoured its wide use in the calculation of expected rate of Dark Matter particle- 
nucleus interaction, it doesn't match with astrophysical observations and presents an 
unphysical behaviour. In fact, the density profile has a singularity in the origin and 
implies a total infinite mass of the halo unless introducing some cut-off at large radii. 
Thus, the use of more realistic halo models is mandatory since the model dependent 
results significantly vary. An extensive discussion about some of the more credited 
realistic halo models has been reported in ref. [17, 18] and have been considered in the 
results given in the following (see Table 2). In particular, the considered halo model 
classes correspond to: i) spherically symmetric matter density with isotropic velocity 
dispersion (A); ii) spherically symmetric matter density with non- isotropic velocity 
dispersion (B); iii) axisymmetric models (C); iv) triaxial models (D); v) moreover, in 
the case of axisymmetric models it is possible to include either an halo co-rotation or an 
halo counter-rotation. The parameters of each halo model have been chosen taking into 
account the available observational data. Thus, considering the allowed range for the 
local velocity of Dark Matter particles vq = (220 ±50) km (90% C.L.), the allowed 
range of local density po has been evaluated [17] taking into account the following 
physical constraints: i) the amount of flatness of the rotational curve of our Galaxy, 
considering conservatively 0.8 -wq ^ v^ot ^ 1-2 -fO) where v^^^ is the value of rotational 
curve at distance of 100 kpc from the galactic center; ii) the maximal non dark halo 
components in the Galaxy, considering conservatively 1 • 10^°Mq < M^s < 6 • 1O"'^''M0 
[44, 45] . Although a large number of self-consistent galactic halo models, in which the 
variation of the velocity distribution function is originated from the change of the halo 
density profile or of the potential, have been considered, still many other possibilities 
exist. 

The proper knowledge of other quantities is also necessary such as e.g. the re- 
coil/electron response ratio for the given nucleus in the given detector and energy 
range, (named quenching factor). Of course, significant differences are often present 
in literature for the measured value of this recoil /electron response ratio even for the 
same nucleus in the same kind of detector as shown in ref. [18]. 

In conclusion, just as a corollary of the model independent result over the seven 
annual cycles, in the following some of the many possible model dependent quests for 
a Dark Matter candidate are summarized. They have been obtained by considering 
the halo models previously mentioned for three of the possible values of the local 
velocity vq: 170 km/s, 220 km/s and 270 km/s and for the halo density values as in 
the prescriptions of ref. [17, 18]. The escape velocity has been maintained at the fixed 
value: 650 km/s. Of course, it is worth to note that the present existing uncertainties 
affecting the knowledge of the escape velocity will significantly extend allowed regions 
e.g. in the cases of preferred inelastic WIMPs and of light mass WIMP candidates; 
its effect would be instead marginal at large WIMP masses. All these scenarios have 
been investigated in some discrc;te cases cither considciring the mean values of the 
parameters of the used nuclear form factors and of the measured quenching factors 
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Table 2: Summary of the considered consistent halo models [17, 18]. The labels in 
the first column identify the models. In the third column the values of the related 
considered parameters are reported [17, 18]; other choices are also possible as well as 
other halo models. The models of the Class C have also been considered including 
possible co-rotation and counter-rotation of the dark halo. 



Class A: spherical pw^ isotropic velocity dispersion 


AO 


Isothermal Sphere 








Al 


Evans' logarithmic 






-Rc = 5 kpc 


A2 


Evans' power-law 






Rc = 16 kpc, 13 = 0.7 


A3 


Evans' power-law 






Rc = 2 kpc, f3 = -0.1 


A4 


Jaffe 






a = 1, ^ = 4, 7 = 2, a = 160 kpc 


A5 


NFW 






a = 1, /? = 3, 7 = 1, a = 20 kpc 


A6 


Moore et al. 






a = 1.5, /3 = 3, 7 = 1.5, a = 28 kpc 


A7 


Kravtsov et al. 






a = 2, /3 = 3, 7 = 0.4, a = 10 kpc 


Class B: spherical pw? no n— isotropic velocity dispersion 


(Osipkov-Merrit, f3o = 0.4) 








Bl 


Evans' logarithmic 






-Rc = 5 kpc 


B2 


Evans' power-law 






Rc = 16 kpc, /? = 0.7 


B3 


Evans' power-law 






-Rc = 2 kpc, f3 = -0.1 


B4 


Jaffe 






a = 1, ^ = 4, 7 = 2, a = 160 kpc 


B5 


NFW 






a = 1, /3 = 3, 7 = 1, a = 20 kpc 


B6 


Moore et al. 






a = 1.5, /3 = 3, 7 = 1.5, a = 28 kpc 


B7 


Kravtsov et al. 






a = 2, /? = 3, 7 = 0.4, a = 10 kpc 


Class C: Axisymmetric pw 


CI 


Evans' logarithmic 






-Rc = 0, g = 1/V2 


C2 


Evans' logarithmic 






-Rc = 5 kpc, q = 1/V2 


C3 


Evans' power-law 






Rc = 16 kpc, q = 0.95, /? = 0.9 


C4 


Evans' power-law 






Rc = 2 kpc, q = I/V2, 13 = -0.1 


Class D: Triaxial pw (q = 0.8, 


P = 


0.9) 


Dl 


Earth on maj. axis, rad. anis. 




5 = -1.78 


D2 


Earth on maj. axis, tang. anis. 




(5 = 16 


D3 


Earth on interm. axis, rad. 


anis. 




(5 = -1.78 


D4 


Earth on interm. axis, tang 


. anis. 


5 = 16 



(case A) or adopting the same procedure as in refs. [15, 16]. The latter one has 
been obtained by varying either: i) the mean values of the measured ^^Na and ^^''I 
quenching factors [32] up to +2 times the errors; ii) the nuclear radius, r„, and the 
nuclear surface thickness parameter, s, in the SI Form Factor [46] from their central 
values down to -20%; iii) the b parameter in the considered SD form factor from the 
given value [47] down to -20% (case B). Moreover, we have also considered one of 
the possible more extreme cases where the Iodine nucleus parameters are fixed at the 
values of case B, while for the Sodium nucleus one considers: i) ^^Na quenching factor 
at the lowest value measured in literature; ii) the nuclear radius, r„, and the nuclear 
surface thickness parameter, s, in the SI Form Factor [46] from their central values up 
to -1-20%; iii) the h parameter in the considered SD form factor from the given value 
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[47] up to +20% (case C). 

Finally, no restriction on the mass of the Dark Matter particle has been adopted in 
these analyses; hence, we have just marked on those figures the lowest bound on the 
neutralino mass as derived from the LEP data in the adopted supersymmetric schemes 
based on GUT assumptions [48]. In fact, other model assumptions are possible and 
would imply significant variations of some accelerators bounds, allowing neutralino 
mass down to 6 GeV (see e.g. the recent refs. [49, 50, 51]); in addition, other low 
mass candidates can be considered as well. It is worth to note that the LEP model 
dependent mass limit - when considered - selects the WIMP-Iodine elastic scatterings 
as dominant because of the adopted scaling laws and of kinematical arguments, while 
DAMA/Nal is intrinsically sensitive both to low and high WIMP mass having both a 
light (the ^^Na) and a heavy (the ^^^I) target-nucleus. 

The results presented by DAMA/Nal on the corollary quests for the candidate 
particle over the seven annual cycles are calculated taking into account the time and 
energy behaviours of the single-hit experimental data. In particular, the likelihood 
function requires the agreement: i) of the expectations for the modulated part of the 
signal with the measured modulated behaviour for each detector and for each energy 
bin; ii) of the expectations for the unmodulated component of the signal with the 
respect to the measured differential energy distribution and - since ref. [13] - also 
with the bound on recoils obtained by pulse shape discrimination from the devoted 
DAMA/NaI-0 data[32]. The latter one acts in the likelihood procedure as an experi- 
mental upper bound on the unmodulated component of the signal and - as a matter 
of fact - as an experimental lower bound on the estimate of the background levels by 
the maximum likelihood procedure. Thus, the C.L.'s, we quote for allowed regions, 
already account for compatibility with the measured differential energy spectrum and 
with the measured upper bound on recoils. In particular, in the following for sim- 
plicity, the results of these corollary quests for the candidate particle are presented in 
terms of allowed regions obtained as superposition of the configurations corresponding 
to likelihood function values distant more than 4fT from the null hypothesis (absence 
of modulation) in each of the several (but still a limited number) of the possible model 
frameworks considered here. Obviously, these results are not exhaustive of the many 
scenarios possible at present level of knowledge (e.g. for some other recent ideas see 
[22, 42, 52]) and larger sensitivities than those reported in the following would be 
reached when including the effect of other existing uncertainties on assumptions and 
related parameters [18]. 

In the most general scenario - to which the DAMA/Nal target nuclei are fully 
sensitive - both the SI and the SD components of the cross section are present. In 
this general scenario the data give an allowed volume in the 4-dimensional space {mw, 
£,(Jsi, S,(JsD, 6)- Fig. 7 just shows slices of this 4-dimensional allowed volume in the 
plane ^asi vs ^asD for some of the possible 6 and mw values for the considered model 
frameworks. We just note that experiments using either even-spin target nuclei (as Ar 
and most of Ge, Xe, Te isotopes) or odd-spin Ge, Xe or Te isotopes cannot explore 
most of the allowed volume. Prom the given figures it is clear that at present either 
a purely SI or a piu-ely SD or mixed SI&SD configurations are compatible with the 
experimental data of the seven annual cycles. 
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Figure 7: Case of a WIMP with mixed SI&SD interaction in the given model frame- 
works. Coloured areas: example of slices (of the 4-dimensional allo'wed volume) in the 
plane £,asi vs S,<jsd for some of the possible mw and 6 values. Four SD couplings are 
reported as examples: i) = (a„ =0 and or \ap\ » |a„|) corresponding to 
a particle "with null SD coupling to neutron; ii) 9 = ^/4 (ap = an) corresponding to a 
particle with the same SD coupling to neutron and proton; iii) 9 = Tr/2 (a„ ^ and 
ttp = or |a„| >> \ap\) corresponding to a particle -with null SD couplings to proton; 
iv) 9 = 2.435 rad = -0.85) corresponding to a particle with SD coupling through 
Zq exchange. The case Op = — a„ is nearly similar to the case iv). Inclusion of other 
existing uncertainties on parameters and models would further extend the regions; 
for example, the use of more favourable form factors and/or of more favourable spin 
factors than the considered ones would move them towards lower cross sections. 



Often the purely SI interaction with ordinary matter is assumed to be dominant 
since e.g. most of the used target-nuclei are practically not sensitive to SD interac- 
tions (instead, ^^Na and ^^^I nuclei are fully sensitive) and the theoretical calculations 
and comparisons are even much more complex and uncertain. Therefore, following 
the analogous procedure as the general case, we have exploited for the same model 
frameworks the purely SI scenario, obtaining the allowed region in the plane m^^ and 
^asi shown in Fig. 8 - left. Of course, best fit values of cross section and WIMP mass 
span over a large range in the considered model frameworks. 

Moreover, configurations with ^cs/ even much lower than those shown in Fig. 8 - 
left are accessible in case an even small SD contribution is present in the interaction 
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Figure 8: On the left : Case of a WIMP with dominant SI interaction for the given 
model frameworks. Region allowed in the plane {mw, ^<ysi)- The vertical dotted line 
represents a bound in case of a neutralino candidate when supersymmetric schemes 
based on GUT assumptions are adopted to analyse the LEP data; the low mass region 
is allowed for neutralino when other schemes are considered and for every other WIMP 
candidate; see text. While the area at WIMP masses above 200 GeV is allowed only 
for few configurations, the lower one is allowed by most configurations (the colored 
region gathers only those above the vertical line) [18]. The inclusion of other existing 
uncertainties on parameters and models would further extend the region; for example, 
the use of more favourable SI form factor for Iodine alone would move it towards lower 
cross sections. On the right : Example of the effect induced by the inclusion of a SD 
component different from zero on allowed regions given in the plane S,asi vs mw ■ In 
this example the Evans' logarithmic axisymmetric C2 halo model with vo — 170 km/s, 
Po equal to the maximum value for this model and a given set of the parameters' 
values (see [18]) have been considered. The different regions refer to different SD 
contributions for the particular case oi 9 — 0: (tsd — (a), 0.02 pb (b), 0.04 pb 
(c), 0.05 pb (d), 0.06 pb (e), 0.08 pb (f). Analogous situation is found for the other 
model frameworks. 

as shown as in an example in Fig. 8 - right. Analogous situation is found for other 
model frameworks. A comparison of the DAMA/Nal purely SI allowed region (given 
in Fig. 8 - left) with the theoretical expectations for a purely SI coupled neutralino 
candidate in a MSSM with gaugino mass unification at GUT scale released is shown 
in Fig. 9 as taken from ref [53] . 

Analogously, one can consider the pure SD coupling. In this scenario one obtain 
an allowed volume in the 3-dimensional space {mw, ^o'sd, 0). Just examples of some 
slices of this allowed volume at given 9 is shown in Fig. 10. Considerations similar to 
the first case hold. 

Finally, also the inelastic Dark Matter particle scenario has been analised obtaining 
an allowed volume in the 3-dimensional space (^dp, mw, ^)- For simplicity. Fig. 11 
shows just few slices of such an allowed volume at some given WIMP masses. We 
remind that in these calculations Vesc has been assumed at fixed value (as in the 
previous cases) , while its present uncertainties can play a significant role in this scenario 
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Figure 9: Figure taken from ref. [53]: theoretical expectations of £,(tsi versus mw in 
the purely SI coupling for the particular case of a neutralino candidate in MSSM with 
gaugino mass unification at GUT scale released; the curve surrounds the DAMA/Nal 
purely SI allowed region as in Fig. 8- left. 
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Figure 10: Case of a WIMP with dominant SD interaction in the given model frame- 
works. Examples of slices (of the 3-dimcnsional allowed volume) in the plane {mw, 
^(Jsd) with 9 = and 9 = 2.435 {Zq coupling). For the definition of the vertical 
line and of the coloured area see previous figure caption. Inclusion of other existing 
uncertainties on parameters and models would further extend the SD allowed regions. 
For example, the use of more favourable SD form factors and/or more favourable spin 
factors would move them towards lower cross sections. Values of S,asD lower than 
those corresponding to these allowed regions are possible also e.g. in case of an even 
small SI contribution. 



of WIMP with preferred inelastic scattering. 

In conclusion, the data of the seven annual cycles have also allowed corollary inves- 
tigations on the features of the Dark Matter particle candidate in some of the many 
possible scenarios. However, although several scenarios have been investigated, these 
corollary analyses are not exhaustive at all because of the present poor knowledge on 
many astrophysical, nuclear and particle Physics needed assumptions. Thus, further 
efforts on the topics are in progress and a relevant contribution is expected by the new 
DAMA/LIBRA set-up, now running. 
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Figure 11: Case of a WIMP with preferred inelastic interaction in the given model 
frameworks. Examples of slices (coloured areas) of the allowed volumes (^Cp, (5, raw) 
for some mw values for WIMP with preferred inelastic interaction. Inclusion of other 
existing uncertainties on parameters and models would further extend the regions; for 
example, the use of a more favourable SI form factor for Iodine and different escape 
velocity would move them towards lower cross sections [18]. 

5 The present situation in the field 

5.1 Direct detection 

As already mentioned, no other experiment, whose result can be directly compared in 
a model independent way with that of DAMA/Nal, is available so far in the field of 
Dark Matter direct detection. 

In fact, most of the activities, started in the 90's, are still at R&D stage and/or have 
released marginal exposures with the respect to the many years of existence and to the 
several used detectors. This is the case of CDMS and EDELWEISS experiments, while 
the Zeplin experiment is more recent [54, 55, 56]. Since these experiments have claimed 
to have "excluded" DAMA/Nal, we will briefly point to the attention of the reader 
only few arguments. In particular. Table 3 summarizes some items for comparison. 

Firstly, let us preliminarily assume as fully correct the " selected" number of events, 
the energy threshold, the energy scale, etc. quoted by those experiments (see Table 3) 
and let us consider if - at least under this hypothesis - their claims might be justified. 
The answer is obviously not; in fact: i) they give a single model dependent result 
using ""'Ge or "°*Xe target, while DAMA/Nal gives a model independent result using 
^■^Na and ^^^I targets; ii) in the single (of the many possible) model scenario, they 
consider, they "fix" all the astrophysical, nuclear and particle physics assumptions at 
a single choice; the same is even for the experimental and theoretical parameters values 
needed in the calculations. In addition, DAMA/Nal is generally quoted there in an 
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Tabic 3: Features of the DAMA/Nal results on the WIMP annual modulation signa- 
ture over the seven annual cycles [18] with those of refs. [54, 55, 56]. 





DAMA/Nal 


CDMS-II 


Edelweiss- 1 


Zeplin-I 


Signature 


Annual 
modulation 


None 


None 


None 


Target-nuclei 




""'Ge 


""'Ge 


""'Xe 


Technique 


well known 


poorly 
experienced 


poorly 
experienced 


critical optical 
liquid/gas interface 
in this realization 


Target mass 


~ 100 kg 


0.75 kg 


0.32 kg 


~ 3 kg 


Exposure 


~ (1.1 ■ lO'^) kg ■ day 


19.4 kg ■ day 


30.5 kg ■ day 


280 kg ■ day 


Depth of the 
experimental site 


1400 m 


780 m 


1700 m 


1100 m 


Software energy 
threshold 


2 keV e.e. 
(5.5 - 7.5 p.e./keV) 


10 keV e.e. 


20 keV e.e. 


2 keV e.e. 
(but: a/E = 100% 
mostly 
1 p.c./keV; [56]) 

(2.5 p.e./keV 
for 16 days; [57]) 


Quenching 
factor 


Measured 


Assumed = 1 


Assumed = 1 
(see also [58]) 


Measured 


Measured event 
rate in low 
energy range 


~ 1 cpd/kg/keV 


??, claimed 7's 
larger than CDMS-I 
{c^ 60 cpd/kg/keV, 
10^ events) 


2± 10^ events 
total 


~ 100 cpd/kg/keV 


Claimed events 
after rejection 
procedures 




either or 1 


2 (claimed taken 
in a noisy period!) 


~ 20-50 cpd/kg/keV 
after rejection and 
?? after standard PSD 
[56, 57] 


Events satisfying 
the signature 
in DAMA/Nal 


modulation 
amplitude 
integrated over the 
given exposure 
~ 10^ events 


insensitive 


insensitive 


insensitive 


Expected number 

of events from 
DAMA/Nal effect 




from few down 
to zero depending 

on the models 
(and on quenching 
factor) 


from few down 
to zero depending 

on the models 
(and on quenching 
factor) 


depends on 
the models 
(even zero) 



uncorrect, partial and unupdated way and existing scenarios to which DAMA/Nal is 
fully sensitive - on the contrary of the others - are ignored. 

Let us now briefly comment also some of the experimental aspects. In particular, 

the counting rate of the Gc bolometers experiments is very high and few/zero events 
are claimed after applying several strong and hardly safe rejection procedures (involv- 
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ing several orders of magnitude). They usually claim to have an "event by event" 
discrimination between noise + electromagnetic background and recoil + recoil-like 
(neutrons, end-range alphas, fission fragments,...) events by comparing the bolometer 
and the ionizing signals for each event, but their results are, actually, largely based 
on " a priori" huge data selections and on the application of other preliminar rejection 
procedures (such as e.g. the one on the so-called surface electrons), which are generally 
poorly described and often not completely quantified. Moreover, most efficiencies and 
physical quantities entering in the interpretation of the claimed selected events have 
never been properly accounted; as an example, we mention the case of the bolome- 
ter quenching factor of the recoil target nuclei. In fact, for the bolometer signals the 
quenching factor (on which the energy threshold and the energy scale rely and, hence, 
also the claimed sensitivity for the given model dependent exclusion plots) is arbitrarily 
assumed to be exactly equal to one. Up to now, only one measurement has been made 
available for a given detector [58]; it offers the value: 0.87 ± 10%{stat.) ± 10%{syst.), 
which is - within the error - compatible with one. but - at the same time - also 
compatible with much smaller values. Thus, any bolometer result, obtained without 
considering e.g. the uncertainties about the unknown value of the quenching factor 
and, hence, about the energy threshold and energy scale, has to be considered partial 
and arbitrary. For completeness we also mention that the reproducibility of the results 
over different running periods has not been proved as well as the values of the effec- 
tive sensitive volumes for the read-outs of the two signals for each event and related 
quantities; obviously, further uncertainties are present when, as done in some cases, a 
neutron background modeling and subtraction is pursued in addition. 

As regards Zeplin-I [56, 57] , a very low energy threshold is claimed (2 kcV) , although 
the light response is very poor: between ~ 1 ph.e./keV [56] (for most of the time) and 
~ 2.5 ph.e./keV (claimed for 16 days) [57] ^. Moreover, a strong data filtering is 
applied to the high level of measured counting rate (see Table 3) by hardware vetoes, 
by fiducial volume cuts and, largely, by applying down to few keV a standard pulse 
shape discrimination procedure, although the LXe scintillation pulse profiles (pulse 
decay time < 30 ns) are quite similar to the PMT noise events in the lower energy 
bins and in spite of the poor light response. Quantitative information on experimental 
quantities related to the used procedures has not yet been given [56, 57]. 

In conclusion, those claims for contradiction have intrinsecally no scientific bases. 

5.2 Indirect detection 

The Dark Matter particles, via their annihilation either in the celestial bodies (such 
as Earth and Sun) or in the Galactic halo could give rise to high energy neutrinos, 
positrons, antiprotons and gamma rays. Therefore, they could be indirectly detected by 
looking either for "upgoing" muons - produced by - in underground, underwater 
or under-ice detectors or for antimatter and gamma rays in the space. However, it 
is worth to remark that no direct model independent comparison can be performed 
between the results obtained in direct and indirect searches. 

^For comparison we remind that the data of the DAMA/LXe set-up, which has a similar light 
response, are analysed by using the much more realistic and safer software energy threshold of 13 keV 
[25]. 
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In the case of "upgoing" muons in terrestrial detectors, the expected n flux is 
the key quantity. However, several sources of uncertainties are present in the related 
estimates (and, therefore, in the obtained results) such as e.g. the assumption that 
a "steady state" has been reached in the considered celestial body and the estimate 
and subtraction of the existing competing processes, offered by the atmospheric neu- 
trinos. Model dependent analyses with a similar approach have been carried out by 
large experiments deep underground such as e.g. MACRO and Super kamiokande. In 
particular, the case of the neutralino candidate in MSSM has been discussed in [59], 
showing that their model dependent results were not in conflict with DAMA/Nal. 

As we mentioned, the annihilation of the Dark Matter particles in the galactic halo 
could also produce antimatter particles and gamma rays. The antimatter searches have 
to be carried out outside the atmosphere, i.e. on balloons or satellites. In particular, 
the Dark Matter particles annihilation would result in an excess of antiprotons or 
of positrons over an estimated background arising from other possible sources. The 
estimate and subtraction of such a background together with the influence of the Earth 
and of the galactic magnetic field on these particles plays a crucial role on the possibility 
of a reliable extraction of a signal. However, at present an excess of positrons with 
energy ~ 5 — 20 GeV has been suggested by [60] and other experiments. Interpreted 
in terms of Dark Matter particles annihilation it gives a result not in conflict with the 
effect observed by DAMA/Nal [60]. 

As regards the possibility to detect 7's from Dark Matter particles annihilation in 
the galactic halo, experiments in space are planned. However, at present it is difficult 
to estimate their possibilities considering e.g. the background level, the uncertainties 
in its reliable estimate and subtraction as well as the smallness of the expected signal 
(even more, if a subdominant component would be present) when properly calculated 
with rescaling procedure. However, in ref. [61, 62] the presence of a 7 excess from 
the center of the Galaxy in the EGRET data [63] has already been suggested. This 
excess match with a possible Dark Matter particles annihilation in the galactic halo 
[61, 62] and is not in conffict with the DAMA/Nal model independent result previously 
published. 

For completeness, we remind that recently it has been suggested [64] that the pos- 
itive hints from the indirect detection - namely the excess of positrons and of gamma 
rays in the space - and the effect observed by DAMA/Nal can also be described in 
a scenario with multi-component Dark Matter in the galactic halo, made of a sub- 
dominant component of heavy neutrinos of the 4*'' family and of a sterile dominant 
component. In particular (see Fig. 12), it has been shown that an heavy neutrino 
with mass around 50 GeV can account for all the observations, while the inclusion 
of possible dumpiness of neutrino density as well as new interactions in the heavy 
neutrino annihilation, etc. can lead to wider mass ranges: from about 46 up to about 
75 GeV (see ref. [64] for details). 
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Figure 12: Figures taken from ref. [64]: Case of a subdominant heavy 4*'* neutrino 
candidate in the plane local density fraction versus the heavy neutrino mass. The 
favorable region obtained from the DAMA/Nal data (grey dashed line when using 
the Evan's halo model; solid line when using the other halo models) and the best- 
fit density parameters deduced from cosmic gamma-radiation (from halo and galactic 
center), positron and antiproton analysis are shown (left panel). The effect of the 
inclusion of possible neutrino dumpiness is also reported (right panel). See ref. [64] 
for details. 

6 Toward the future: from DAMA/Nal to 
DAMA/LIBRA and beyond 

The large merits of highly radiopure Nal(Tl) set-up have been demonstrated in the 
practice by DAMA/Nal which has been the highest radiopure set-up available in the 
field, has effectively pursued a model independent approach to investigate Dark Matter 
particles in the galactic halo collecting an exposure several orders of magnitude larger 
than those available in the field and has obtained many other complementary or by- 
products results. 

In 1996 DAMA proposed to realize a ton set-up [65] and a new R&D project for 
highly radiopure Nal(Tl) detectors was funded at that time and carried out for several 
years in order to realize as an intermediate step the second generation experiment, 
successor of DAMA/Nal, with an exposed mass of about 250 kg. 

Thus, new powders and other materials have been selected, new chemical/physical 
radiopurification procedures in Nal and Til powders have been exploited, new grow- 
ing/handling protocols have been developed and new prototypes have been built and 
tested. As a consequence of the results of this second generation R&D, the new ex- 
perimental set-up DAMA/LIBRA (Large sodium Iodide Bulk for RAre processes), 
~250 kg highly radiopure Nal(Tl) crystal scintillators (matrix of twenty-five ~ 9.70 kg 
Nal(Tl) crystals), was funded and realised. After the completion of the DAMA/Nal 
data taking in July 2002, the dismounting of DAMA/Nal occurred and the installa- 
tion of DAMA/LIBRA started. In particular, the experimental site as well as many 
components of the installation itself have been implemented (environment, shield of 
PMTs, wiring, HP Nitrogen system, cooling water of air conditioner, electronics and 
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Figure 13: Left: the installation of the 25 Nal(Tl) crystals (9.70 kg each one) of 
DAMA/LIBRA in HP Nitrogen atmosphere. Right: One of the final stages of the 
detectors' installation. Bottom: Picture of the photomultiplier shield made of highly 
radio-pure Copper. All the used materials have been deeply selected for radiopurity 
(see for example the cables with teflon envelop). All the procedures as well as these 
photos have been carried out in HP Nitrogen atmosphere. 

DAQ, etc.). In particular, all the Cu parts have been chemically etched before their 
installation following a new devoted protocol and maintained in HP Nitrogen atmo- 
sphere until the installation. All the procedures performed during the dismounting of 
DAMA/Nal and the installation of DAMA/LIBRA detectors have been carried out in 
HP Nitrogen atmosphere (see Fig. 13). 

DAMA/LIBRA is taking data since March 2003 and the first data release will, 
most probably, occur when an exposure larger than that of DAMA/Nal will have been 
collected and analysed in all the aspects. Just as an example of the quality of the 
data taking, Fig. 14 shows the stability of the calibration factor and of the ratio of 
the peaks' positions of the ^'''^Am source during about one year of data taking. 

The highly radiopure DAMA/LIBRA set-up is a powerful tool for further inves- 
tigation on the Dark Matter particle component in the galactic halo having all the 
intrinsic merits already mentioned in section 2 and a larger exposed mass, an higher 
overall radiopurity and improved performances with the respect to DAMA/Nal. Thus, 
DAMA/LIBRA will further investigate the 6.3 a C.L. model independent evidence 
pointed out by DAMA/Nal with increased sensitivity in order to reach even higher 
C.L.. Moreover, it will also offer an increased sensitivity to improve corollary quests 
on the nature of the candidate particle, trying to disentangle at least among some of 
the many different possible astrophysical, nuclear and particle physics models as well 
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Figure 14: An example of the stability of the calibration factor (tdcal) and of the ratio 
of the peaks' positions (a) of the measured energy distribution of the ^'^^Am source 
during about one year of data taking. 



as to investigate other new possible scenarios. 

In the following some of the main topics - not yet well known at present and which 
can affect whatever model dependent result and comparison - will be mentioned. They 
will be addressed by the highly radiopure DAMA/LIBRA in a way often unique and 
always reliable, cheap and competitive. They are: 

• the velocity and spatial distribution of the Dark Matter particles in the galac- 
tic halo. It has been shown that the naive description of the galactic halo as 
an isothermal halo is an unphysical and non-realistic approximation which sig- 
nificantly affects model dependent evaluations (exclusion plots, allowed regions, 
etc.) and comparisons. Other modelings (not exhaustive at all), many of them 
based on N-bodies simulations, have been considered in literature and some of 
them have been discussed at some extent in [17, 18] and references therein. Some 
of these models could be significantly discriminated by DAMA/LIBRA. 

• the effects induced on the Dark Matter particles distribution in the galactic 
halo by contributions from satellite galaxies tidal streams. Recently it has been 
pointed out [22] that contributions to the Dark Matter particles in the galactic 
halo should be expected from tidal streams from the Sagittarius Dwarf elliptical 
galaxy. Considering that this galaxy was undiscovered until 1994 and consider- 
ing galaxy formation theories, one has to expect that also other satellite galaxies 
do exist and contribute as well. In particular, the Canis Major satellite Galaxy 
has been pointed out as reported in 2003 in ref. [66]; it can, in principle, play a 
very significant role being close to our galactic plane. At present, the best way 
to investigate the presence of a stream contribution is to determine in accurate 
way the phase of the annual modulation, to, as a function of the energy; in fact. 
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for a given halo model to would be expected to be (slightly) different from 152.5 
d and to vary with energy (see Fig. 15). 
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Figure 15: Expected behaviours of the phase, to, of the annual modulation signal as 
function of the energy when considering: i) only galactic halo ("no Sgr"); ii) galactic 
halo (C2 halo model with vq = 220 km/s, po equal to the maximum value for this 
model) and a contribution from Sagittarius Dwarf galaxy ("C2"); iii) galactic halo 
(A5 halo model with vq = 220 km/s, po equal to the maximum value for this model) 
and a contribution from Sagittarius Dwarf galaxy ("A5"). The contributions from 
Sagittarius Dwarf galaxy have been taken in both cases with a density equal to 4% of 
Pq. The light shadow region is the final result of DAMA/Nal on the to value for the 
cumulative energy interval (2 - 6) keV, while the dark shadow region is the expectation 
on to assuming an experiment with the same features as DAMA/Nal, an exposure of 
3 • 10^ kg • day and the same central value for to- 

• the effects induced on the Dark Matter particles distribution in the galactic halo 
by the existence of caustics. It has been shown that the continuous infall of Dark 
Matter particles in the galactic gravitational field can form caustic surfaces and 
discrete streams in the Dark Matter particles halo [52]. The phenomenology 
to point out a similar scenario is analogous to that in the previous item; thus, 
DAM A/LIBRA can as well test this possibility. 

• the detection of possible "solar wakes". As an additional verification of the pos- 
sible presence of contributions from streams of Dark Matter particles in our 
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galactic halo, DAMA/LIBRA can investigate also the gravitational focusing ef- 
fect of the Sun on the Dark Matter particle of a stream. In fact, one should 
expect two kinds of enhancements in the Dark Matter particles flow: one named 
"spike", which gives an enhancement of Dark Matter particle density along a 
line collinear with the direction of the incoming stream and of the Sun, and 
another, named "skirt", which gives a larger Dark Matter particle density on 
a surface of cone whose opening angle depends on the stream velocity. Thus, 
DAMA/LIBRA will investigate such a possibility with high sensitivity through 
second-order time-energy correlations. 

• the study of possible structures as dumpiness with small scale size. Possible 
structures as dumpiness with small scale size could, in principle, be pointed out 
by exploiting a large exposure which can be collected by DAMA/LIBRA. 

• the coupling(s) of the Dark Matter particle with the ^^Na and ^^^/ and its nature. 
As mentioned, several large uncertainties are linked to the coupling(s) between 
the Dark Matter particle and the target-nuclei. DAMA/LIBRA, exploiting a new 
large exposure, will allow to better constrain the related aspects. In addition, 
analyses in model frameworks suitable for e.g. mirror Dark Matter [67] and for 
particles from multi-dimensional Kahiza-Klcin like theories will be performed as 
well and compared with the other scenarios analised so far. 

• scaling laws and, cross sections. At present just simple scaling laws are used to 
scale all the nuclear cross sections to a common nucleon cross section; however, 
they are a large source of uncertainties in model dependent results and com- 
parisons. For example, recently, it has been pointed out [42] that, even for the 
neutralino candidate, these assumptions (which hold in the case of model with 
one-nucleon current) are arbitrary when two-nucleon current with pion exchange 
are introduced. Thus, the presence of two target-nuclei in the Nal(Tl) detectors 
of DAMA/LIBRA could in principle offer a probe for that in the large exposure 
this experiment will collect. 

A large work will be faced by DAMA/LIBRA, which is in addition the intrinse- 
cally most sensitive experiment in the field of Dark Matter because of its radiopurity. 
exposed mass and high duty cycle. These qualities will also allow DAMA/LIBRA to 
further investigate with higher sensitivity several other rare processes. 

Finally, at present a third generation R&D effort toward the possible Nal(Tl) ton 
set-up has been funded and related works have already been started. 

7 Conclusion 

DAMA/Nal has been a pioneer experiment running at LNGS for several years and 
investigating as first the model independent WIMP annual modulation signature with 
suitable sensitivity and control of the running parameters. During seven independent 
experiments of one year each one, it has pointed out at 6.3 a C.L. the presence of a 
modulation satisfying all the many peculiarities of a WIMP induced effect. Neither 
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systematic effects nor side reactions able to account for the observed modulation am- 
plitude and to contemporaneously satify all the requirements of the signature have 
been found. DAMA/Nal has also pointed out the complexity of corollary investiga- 
tions on the nature of the candidate particle, because of the present poor knowledge 
on the many astrophysical, nuclear and particle physics aspects. 

A second generation experiment DAMA/LIBRA has been realized and put in oper- 
ation since March 2003. This new set-up, having a larger exposed mass and an higher 
overall radiopurity, will offer a significantly increased sensitivity to further investigate 
the Dark Matter particle component in the galactic halo pointed out at 6.3 a C.L. 
by DAMA/Nal. Moreover, further efforts towards the creation of ultimate radiopure 
Nal(Tl) set-ups are already in progress; in particular, a third generation R&D is in 
progress towards the possible 1 ton set-up we proposed in 1996. 

References 

[1] F. Zwicky, Helv. Phys. Acta 6 (1933) 110 
[2] S. Smith, Astrophys. J. 83 (1936) 23 

[3] V.C. Rubin and W.K. Ford, Astrophys. J. 159 (1970) 379; M. Roberts and A.H. 

Rots, Astron. Astrophys. 26 (1973) 483 
[4] D. N, Spergel et al., astro-ph/ 0302209 

[5] A. Riess et al., Astronom. J. 116 (1998) 1009; E.D. Perlmutter et al., Astrophys. 

J. 517 (1999) 565 
[6] C. L. Bennett et al., astro-ph/0302207 

[7] R.A.C. Croft et al., Phys. Rev. Lett. 83 (1999) 1092; O. Elgaroy et al. astro- 

ph/0204152 

[8] R. Bernabei et al., Phys. Lett. B 515 (2001) 6 

[9] P. Belli, R. Bernabei, C. Bacci, A. Incicchitti, R. Marcovaldi, D. Prosperi, DAMA 
proposal to INFN Scientific Committee II, April 24*'' 1990. 

[10] R. Bernabei et al., Phys. Lett. B 424 (1998) 195 
[11] R. Bernabei et al., Phys. Lett. B 450 (1999) 448 
[12] P. BeUi et al., Phys. Rev. D 61 (2000) 023512 
[13] R. Bernabei et al., Phys. Lett. B 480 (2000) 23 
[14] R. Bernabei et al., Eur. Phys. J. C18 (2000) 283 
[15] R. Bernabei et al., Phys. Lett. B 509 (2001) 197 
[16] R. Bernabei el al., Eur. Phys. J. C23 (2002) 61 
[17] P. BelU et al., Phys. Rev. D 66 (2002) 043503 

[18] R. Bernabei el al.. La Rivista del Nuovo Cimento 26 n.l (2003) 1-73 [astro- 

ph/0307403] 

[19] K.A. Drukier et al., Phys. Rev. D 33 (1986) 3495; K. Freese et al., Phys. Rev. D 
37 (1988) 3388 

[20] P.J.T. Leonard and S. Tremaine, Astrophys. J. 353 (1990) 486; C.S. Kochanek, 
Astrophys. J. 457 (1996) 228; K.M. Cudworth, Astron. J. 99 (1990) 590 



30 



[21] D. Smith and N. Weiner, Phys. Rev. D 64 (2001) 043502 and hep-ph/0402065. 
[22] K.Preese et. al. astro-ph/0309279; Phys. Rev. Lett. 92 (2004) 11301. 
[23] R. Bernabei et al., // Nuovo Cim. A112 (1999) 545 

[24] P. Belli et al., Astrop. Phys. 5 (1996) 217 ; P. Belli et al., II Nuovo Cimento C 19 
(1996) 537 ; P. Belli et al., Phys. Lett. B 387 (1996) 222; Phys. Lett. B 389 (1996) 
783 (err.); P. Belli et al., Phys. Lett. B 465 (1999) 315 ; P. Belli et al., Phys. Rev. 
D 61 (2000) 117301 ; R. Bernabei et al.. New Journal of Physics 2 (2000) 15.1 ; 
R. Bernabei et al., Phys. Lett. B 493 (2000) 12 ; R. Bernabei et al., Nucl. Instr. & 
Methods A 482 (2002) 728 ; R. Bernabei et al., Eur. Phys. J. direct Cll (2001) 1; 
R. Bernabei et al., Phys. Lett. B 527 (2002) 182 ; R. Bernabei et al., Phys. Lett. 
B 546 (2002) 23 . 

[25] R. Bernabei et al., Phys. Lett. B 436 (1998) 379 

[26] R. Bernabei ct al., Astrop. Phys. 7 (1997) 73 ; R. Bernabei et al., II Nuovo Cimento 
A 110 (1997) 189 ; P. Belli et al., Astrop. Phys. 10 (1999) 115; P. Belli et al., 
Nucl. Phys. B 563 (1999) 97 ; R. Bernabei et al., Nucl. Phys. A 705 (2002) 29 ; P. 
Belli et al Nucl. Instr. & Methods A 498 (2003) 352 ; R. Cerulli et al Nucl. Instr. 

& Methods A 525 (2004) 535 . 

[27] M. Ambrosio et al., Astrop. Phys. 7 (1997) 109 
[28] P. Belli et al., II Nuovo Cimento A 101 (1989) 959 
[29] M. Cribicr ct al., Astrop. Phys. 4 (1995) 23 
[30] C. Arpesella et al.. Health Phys. 72 (1997) 629 
[31] M. Wojcik: Nucl. Instrum. & Methods B61 (1991) 8 
[32] R. Bernabei et al., Phys. Lett. B 389 (1996) 757 
[33] R. Bernabei et al., // Nuovo Cimento A 112 (1999) 1541 
[34] P. Belli et al., Phys. Lett. B460 (1999) 236 
[35] R. Bernabei et al., Phys. Lett. B 408 (1997) 439 . 
[36] P. Belli et al., Phys. Rev. C60 (1999) 065501. 
[37] R. Bernabei et al., Phys. Rev. Lett. 83 (1999) 4918 . 
[38] F. Cappclla ct al., Eur. Phys. J.-dircct C14 (2002) 1. 
[39] R. Bernabei et al., Eur. Phys. J. A 23 (2005) 7. 
[40] R. Bernabei et al., ROM2F/2004/30, submitted for publication. 
[41] W.H. Press and G. B. Rybicki, Astrophys. J. 338 (1989) 277; J.D. Scargle, Astro- 
phys. J. 263 (1982) 835 

[42] G. Prezeau et al., Phys. Rev. Lett. 91 (2003) 231301. 
[43] A. Bottino et al., Phys. Lett. B402 (1997) 113 

[44] W. Dehnen and J. Binney, Mon. Not. R. Astron. Soc. 294 (1998) 429 
[45] E. I. Gates, G. Gyuk, M. S. Turner, Phys. Rev. D53 (1996) 4138 
[46] R.H. Helm, Phys. Rev. 104 (1956) 1466 ; A. Bottino et al., Astrop. Phys. 2 (1994) 
77 

[47] M.T. Ressell et al., Phys. Rev. C 56 (1997) 535 
[48] D.E. Groom ct al., Eur. Phys. J. C 15 (2000) 1 

[49] A. Bottino et al., Phys. Rev. D 67 (2003) 063519 ; A. Bottino et al, hep-ph/0304080 



31 



[50] D. Hooper and T. Plehn, MADPH-02-1308, CERN-TH/2002-29, hep-ph/0212226 
[51] G. Belanger, F. Boudjema, A. Pukhov and S. Rosier-Lees, hep-ph/0212227 
[52] F.S. Ling, P. Sikivie, S. Wick, astro-ph/0405231. 

[53] A. Bottino ct al., Phys. Rev. D 69 (2004) 037302 . 

[54] CDMS collaboration, astro-pli/0405033; Phys. Rev. Lett. 84 (2000) 5699 
[55] EDELWEISS collaboration, in the Proc. of NDM03, Japan (2003); Phys. Lett. B 
513 (2001) 15 

[56] N. Smith, talk given at IDM02, York, September 2002 

[57] R. Luscher , talk given at Moriond, march 2003 

[58] E. Simon et al., Nucl. lustrum. & Meth. A507 (2003) 643. 

[59] A. Bottino et al., Phys. Lett. B 402 (1997) 113 ; Phys. Lett. B 423 (1998) 109 ; 
Phys. Rev. D 59 (1999) 095004 iPhys. Rev. D 59 (1999) 095003 ; Astrop. Phys. 10 
(1999) 203 ; Astrop. Phys. 13 (2000) 215 ; Phys. Rev. D 62 (2000) 056006 ; Phys. 
Rev. D 63 (2001) 125003 ; Nucl. Phys. B 608 (2001) 461 

[60] G.L. Kane et al., hep-ph/0108138 

[61] A. MorseUi, talk given at Vulcano 2002 

[62] A. Morselh et al., astro-ph/0211327 

[63] A. Strong et al., Astrophys. J. 537 (2000) 763 

[64] K. Belotsky, D. Fargion, M. Khlopov and R.V. Konoplich, hep-ph/0411093 

[65] R. Bcrnabci ct al., Astrop. Phys. 4 (1995) 45; R. Bernabei, in The identification 

of Dark Matter, World Sc. Pub. (1997) 574. 
[66] R.A. Ibata et al., Man. Not. Roy. Astr. Soc. 348 (2004) 12. 
[67] R. Foot, hep-ph/0308254. 



32 



